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The hnear integral equations defining the Navier-Stokes (NS) transport coefficients for polydis- 
perse granular mixtures of smooth inelastic hard disks or spheres are solved by using the leading 
terms in a Sonine polynomial expansion. Explicit expressions for all the NS transport coefficients are 
given in terms of the sizes, masses, compositions, density and restitution coefficients. In addition, 
the cooling rate is also evaluated to first order in the gradients. The results hold for arbitrary degree 
of inelasticity and are not limited to specific values of the parameters of the mixture. Finally, a 
detailed comparison between the derivation of the current theory and previous theories for mixtures 
is made, with attention paid to the implication of the various treatments employed to date. 



I. INTRODUCTION 

In the first portion II] of this two-paper effort on the development of a kinetic-theory-based description of mixtures, a 
rigorous Chapman-Enskog (CE) expansion of the revised Enskog theory (RET) for an s -component mixture of inelastic 
hard spheres was carried out to first (Navier-Stokes) order in spatial gradients. The expansion was performed about 
a homogeneous cooling state (HCS) , which is the zeroth-order solution of the kinetic equation for the single-particle 
velocity distribution function fi of species i. Unlike previous theories for mixtures which derive from expansions 
about an elastic base state [1, |3, 0, S H, S B S this CE expansion does not impose any constraints on the 
level of dissipation. The resulting theory is thus expected to be applicable to a wide range of restitution coefficients. 
Furthermore, because the derivation used the RET as its starting point (as opposed to the Boltzmann equation as 
used in Ref. [HI), the results are expected to be applicable to dilute and moderately dense systems. The theory is 
not expected to be applicable to systems dense enough for ring collisions to play a significant role, since such velocity 
correlations are not accounted for in the RET. 

This formally exact analysis for Navier-Stokes order hydrodynamics, reported in the companion paper, resulted in 
integral-differential equations for the zeroth-order ff'^ and first-order fl'^^ distribution functions as well as integral 

expressions (in terms of ff^ and fl^^) for each of the equations of state (cooling rate and pressure) and the transport 
coefficients {Dij,Df , Dfj,r], k, A, Dq^ij,Lij}. Of these, only the pressure could be directly evaluated and cast in alge- 
braic expressions of the macroscopic (hydrodynamic) variables. Hence, in this second part of the work, approximate 
methods are used to obtain algebraic equations for the kinetic and coUisional contributions to the cooling rate and 
transport coefficients. In particular, the equations for f^^"^ have previously been solved via a combination of scal- 
ing arguments and an approximation to the distribution function based on leading-order Sonine polynomials. Recent 
results derived for binary granular mixtures at low-density [ill . [l3| have shown that the influence of the non-Gaussian 
(higher-order) corrections of /■^'^■' to the transport coefficients is in general negligible, except for quite large values of 
dissipation and velocities. For this reason, we use leading order in the Sonine polynomial expansion (Maxwellians) at 

different temperatures to evaluate integrals over the distributions fj^°\ The solutions for the f-^^ are also found using 

a truncated polynomial expansion. These forms for f^^"^ and fj^^'^ are then used to obtain practical expressions for 
the cooling rate and transport coefficients. Namely, the resulting, algebraic constitutive equations depend explicitly 
on the hydrodynamic variables only and not on the distribution functions of the mixture. Collision integrals are 
reduced to Gaussian forms and evaluated using standard integration techniques; a computer package of symbolic 
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calculations (Mathematica) was used to check their analytical evaluation. Wherever possible, limiting values of the 
equations of state and transport coefficients were verified by comparison with previously published works in the cases 
of mechanically equivalent particles fll. [l5l | , binary mixtures at low-density [ill, EBl , and simple shear flow states for 
binary mixtures at moderate density \l7\ . All analytic calculations of both papers were performed independently and 
cross-checked. 

Finally, as mentioned in the Introduction of the companion paper, a myriad of treatments has appeared in the 
literature for the derivation of kinetic-theory-based descriptions of granular mixtures. To help put the current effort 
in the context of previous contributions, Sec. VI contains a breakdown of the theoretical contributions to date, along 
with a critical discussion of the ramifications of each treatment. 



II. EULER ORDER PARAMETERS 



The hydrodynamic equations to first order in the spatial gradients (Euler order) have as unknown (phcnomenolog- 
ical) parameters the pressure and the cooling rate, 

p^p{T,{n,}), C(r,{n,}) = C^")(r,{n,;})-fCc/(r,{nJ)V-U. (2.1) 

The pressure p{T, {ui}) has been evaluated exactly in the previous paper [l| [Eqs. (7.18) and (7.19)] with the result 

d/2 * 

p(r, {n,}) =nT+ ^-j- J2 E t^j^^xf^^n.n.il + a,,)T,, (2.2) 

V 2 / i—l j—l 

where henceforth, for the sake of convenience, we will use the notation T, = t/°^ with 7;^°^ defined by Eq. (7.5) of Ref. 
[ij. Detailed forms for the functions Xij^ are discussed in Appendix [C] The leading order cooling rate, C^°H^j {'^iDj 
also is given there by Eq. (7.2) as an integral over f^^"^ . For reasons just mentioned above, these integrals are performed 
using the leading order Sonine approximation for fj^^^ 

/r(v,^/,.,(v,^„.(^)"'=p(-^), (.3) 

and V = V — U is the peculiar velocity. The result is [13, [3] 



C(°)(r,{n,}) = (T,{na) 



^(1 



(2.4) 



where n = ^ ^ is the total density, vq (t) = •\/2T/m is a thermal velocity, m = (^^ rrij )/s, ^ij = rrii/ {rrii + raj ) , and 

9i = niiT/mTi. The C.I'^\t, {rii}) are species cooling rates, measuring the decrease of kinetic energy of each species. 
This is an implicit definition of (T, {ui}). The equality of cooling rates for all i gives s — 1 equations for the species 
temperatures in terms of T. A final equation is given by the condition that the total kinetic energy is the sum of 
species energies nT — X^i^i"-*^!- With the species temperatures determined as functions of T and {n^}, Eq. (|2.4p 
gives the cooling rate C'-^H'^i {"-i})- all of the following expressions, it is understood that the Ti{T, {ui}) have been 
determined in this way. 

At first order in gradients, there is a contribution to the cooling rate from V • U. The proportionality coefficient 
Cc/ is a new transport coefficient for granular fluids. Two different contributions can be identified 

C^ = C^M)+^(i,i). (2.5) 

The coefficient (^^^'^'^ is given by an integral over f^^"^ and its velocity derivative. It has been evaluated explicitly in 
Appendix F of Ref. |1|] with the result 

c^^^°^ = -;|^EE-^".M.^4x!r(i-4)^. (2.6) 
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The second term of Eq. (|2.5|) . C''^'^\ is given by 



c 



(1,1) 



1 T^(d~^)l'i ' 



i±3yEE4^'4'™^M..(l-4) / '^^i / dv2g3/f)(Vi)I?,(V2), (2.7) 



where the unknown functions T)i{W) are the solutions to the hnear integral equations 



Here, Vv = d/dY and the hnear operator C is 



(/:x), = ic'")Vv(vx,) + (Ln, 



where L is the linearized Enskog collision operator 

s 

j=i 



(0) 



7(0) 



(2.8) 



(2.9) 



(2.10) 



J^) [Vi I X„ /,] ^ xf (ri, i)}) 1 dV2 1 dae(5^ • gi2)(^ • gl2) 

X Kr.2x,(V'/)/,(V^') - X.(Vi)/,(V2)] , (2.11) 

and v",V2 are the precoUision velocities associated with vi,V2 for the coUiding pair of species i and j. The same 
hnear operator £ appears below in the definitions of all other transport coefficients as well. The inhomogeneity of 
Eq. jMl), A(T,{nJ), is given by 



D,. 



(0) 



j=i 



drii 



(0)i 



(2.12) 



j=i 



and the operator /Cy./j is given by (|Bip of Appendix |B] An approximate solution to this integral equation is obtained 

by using ff^\ p, ^("^ and Cp''^^ as determined above, and the leading order term in an expansion of I?i(V) in a 
complete set of Sonine polynomials 



I?.(V)^e,^Z5AM(V)F,(V). 
The lowest order Sonine polynomial contributing in this case is 



f.(v, = |||l v-~-v- + 



d + 2m.i 2 d{d + 2) 



(2.13) 



(2.14) 



as determined by the fact that I'i(V) is a scalar and orthogonal to 1 and . Finally, o is the projection of Vi 
along Fi 



ej n = 



d{d + 2) Ui 



dv2?,(V)F,(V). 



(2.15) 



Note that the polynomials are defined for each species in terms of the weight factor /i,M(V) characterized by the 
mass and temperature of that species [ll|. The coefficients ei^n are determined by substituting (|2.13p into the integral 
equation (|2.8[) , multiplying by Fi (V) , and integrating over V (this assures that the integral equation is exactly satisfied 
in the subspace spanned Fi) to get the algebraic equations 



3 = 1 ^ ^ 



(2.16) 
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The collision frequencies il'ij are defined by 

2 



d{d + 2) n. 



J- ( E / rfvF,j(°)[/,MF.,/f ] + I dv^^,ji°)[/f\ ^Mi^,] j , (2.17) 



L [dvF,4\f^'\UMF,], {^^J). (2.18) 



d{d + 2) ni 

The inhomogeneity on the right-hand side of p.l6p is computed in Appendix [B] with the result 



e^,D = / dv F,(y)JC,j,^[dvJ, 



(Oh 
3 I 



= -^^jfjd^ E ^^^Jxif'^tj^^Ji{^ + a^j) [8(d + - 1) 

+4(13 + 2d + 9ay )mj, - ^^ii%6j^{e, + 6'j)(l + ay)^ 

+ 15/ij',07'(e. + 0,f{l + a,,f] . (2.19) 

The explicit form of the collision frequencies ■0^ are displayed in Appendix[Xl so the coefficients Ci^o can be determined 
by solving the algebraic Eqs. (|2.16p with all coefficients determined. 

Finally, is given in terms of Ci^jj by substitution of Eq. (|2.13p into (|2.7p . with the result 



III. NAVIER-STOKES ORDER TRANSPORT COEFFICIENTS 



The hydrodynamic equations to second order in the spatial gradients have contributions to the mass flux, pressure 
tensor, and heat flux from the Chapman-Enskog expansion proportional to gradients of the species densities, tem- 
perature, and flow velocity and to an external applied force (of the same order of magnitude as the gradients). For 
fluid symmetry, this leads to (s — l)(2s -I- 1) transport coefficients for the mass flux, the shear and bulk viscosity 
for the pressure tensor, and 2s^ -I- 1 coefficients for the heat flux. The mass flux coefficients are determined from 
integral equations in the same way as C*'^'^''; just described. The pressure tensor and heat fluxes have "kinetic" and 
"collisional" transfer parts, i.e. of the forms 

The contributions uj^ are determined from integral equations like that for C*"^'^^; the collisional transfer contributions 
have been reduced in the previous paper to specific integrals over and can be computed explicitly by using the 
multi-temperature Maxellians (|2.3p . The analysis is similar to that for the cooling rate and many of the details are 
transferred to the Appendices. 



A. Mass Flux 

To first order in the spatial gradients, the mass flux Jq^^ is given by 

s 

;(i) 



j^:' =-E^^^A,^lnn,-pA^VlnT-^A^F,, (3.1) 



where p — pi — mini is the total mass density of the mixture, Dij are the mutual diffusion coefficients, 
are the thermal diffusion coefRcients and Df^ are the mobility coefficients. These transport coefficients are defined as 

^/dvV.Si(V), (3.2) 



dm A 
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(3.3) 



Df,-~^ld.V-Sl{V). (3.4) 
The functions Ai (V) , Bj (V) , and (V) verify the Hnear integral equations 

^-^C'"^)^) =A„ (3.5) 

{CB^),-n,^A.=Bi, (3.6) 

C + C^°A£') =El, (3.7) 



where the hnear operator C and the inhomogeneous terms A^, B^, and are given by Eqs. (6.16), (6.17), and (6.20) 
of the preceding paper [J] . The approximate solutions to these equations are obtained in the same manner as that 
for ( 12. 8p . The leading order Sonine polynomial in all of these cases is V, since all are vectors 

A(V)^-^i?f/,M(V)V, (3.8) 

m(y) - -^^^^D,,f,,M(y)Y, (3.9) 

£l{^)^-^Dfj,,M{Y)\, (3.10) 

The coefficients are now projections of S^, and £l along V, which are identified in terms of the transport 
coefficients through Eqs. (|3.2p - (|3.4p . Multiplication of Eqs. (|3.5p - (|3.7p by rriiV and integrating over the velocity 
yields the algebraic equations determining all mass flux transport coefficients 

l^^^S^'^m.-J^f^Dj I d.rn.Y.Bliyi (3.12) 



E 

£=1 



I'll 



The new collision frequencies Vij are 



E (^^^ + ^C(°^'5^^) = -\\ d^^i^ ■ Ei(V). (3.13) 



J dvrri^V ■ 4^\f,,Myjl'\ (3.14) 



(3.15) 



Note that the self-collision terms of va arising from [/i,AfV, do not occur in p.l4p since these conserve 
momentum for species i. The above collision frequencies were already evaluated in the Boltzmann limit (except for 
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the factors [HI, [3- 

The details will not repeated here and only the results are quoted in Appendix [X] The 
integrals on the right hand side of Eqs. p. lip . (|3.12p . and p.l3p can be performed using the definitions of A^, B^, 
and . The result is 



dv rriiV ■ A,; = dp— [ 1 — ^ * * 



PiP 



\Y.j dvm,K,/C,,,,[Vv(V/f )], 



(3.16) 



dv TO,- V • B;? 



-dn 



dp 



— [niTi] + d—rij — 
ouj p on 



(3.17) 



dvTOjV • E;? 



(3.18) 



The integrals appearing in Eqs. (j3.16p and (|3.17p that involve the operator /C.y,^[X] have been evaluated in Appendix 
m They are given by Eqs. (|XT8)) . (|B6|) . and With these results, Eqs. ((XTT1) . ((3T2)) . and ((3T3)) become 



puiTi 



(3.19) 



5Z I '^'^ 



mi 



D 



P 



D 



rairrij drij 

df 

1 n£ 



duj rrij drij 



rriimj 



5u + 



'^i-^ — In X, 



(0) 



mi drij 



\nji 



(3.20) 



n^TOi 



TO," 



71,TO,- 



(3.21) 



where ji = Ti/T is the temperature ratio. An explicit form for luj is chosen in Appendix [Cl The solution of the 
set of algebraic equations p.l9p - p.2ip gives the dependence of the coefficients , Dj , and Df^ on the restitution 
coefficients and the composition, the density, and the sizes and masses of the constituents of the mixture. 



B. Pressure Tensor 

The constitutive equation for the pressure tensor P^^, proportional to the velocity gradients, is 

= -77 [dr^Up + drpU^ - ^(5„^V • - kJ^^V • u. (3.22) 

Here, rj is the shear viscosity coefficient and k is the bulk viscosity. The coefficient rj has kinetic and coUisional 
contributions while k only has a coUisional contribution (and so, vanishes for dilute gases) 

r^^n'^+r^", k^k". (3.23) 

The coUisional transfer contributions have been analyzed in the preceding pap er [H [Eqs. (7.21) and (7.22)]. These 
expressions reduce to those previously derived in the monodisperse case [14i Il5| , and in the case of binary mixtures of 
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hard spheres p7| . The integrals over /j^"-' (V) are easily performed in the multi-temperature Maxwellian approximation 
with the results 

E E ;r^-^%-o4+^x!? (1 + a.. ) ^ , (3.24) 



= E E %-4x!°V.(l + + (3.25) 

1. Kinetic contribution rj'' 

As noted above, there is no kinetic part of the bulk viscosity, k'^ = 0, so k is given entirely by p.24p . The kinetic 
contribution to the shear viscosity, ry*^ , is defined by 



^7' = E = - (d-l)id + 2) E / mVaVpC.MV). 



(3.26) 



The second equality identifies the partial contribution rj'^ of the species i to the shear viscosity rj'' in terms of Ci_Q,/3(V), 
which is the solution to the integral equation 



(3.27) 



where Ci^ap is given by Eq. (6.18) in the preceding paper pQ. It is symmetric and traceless so the leading Sonine 
approximation for the function Ci.Q/3(V) is 

C,,a/3(V) ^ -/,.^,(V) Ai?^ (V), (3.28) 

where 

R^,apm = ^^ (VaVp - ^VH^f,^ . (3.29) 

The partial contributions rj^ are obtained by multiplying the integral equation (I3.27P with Ri.ap and integrating over 
the velocity to get the set of equations: 

E(n. - \c^°^S,,)r^' - - E / R^.cpJC^.Advjf], (3.30) 

j=i ^ ^ j=i 

The collision frequencies tu are 



1 



{d- l)(d + 2) n.. 



^ J d^^R^,a^J^?[kMR^,cpJf] + J dv i?,,,^ /,,Mi?^,a/3] j , (3.31) 



1 



[d- l)(d+2)n 



-i;^ /dvi?,,„/3jf [/f ,/,-M%o/3], (i^j). (3.32) 



The explicit forms of tu and are displayed in Appendix \^ and the integral appearing on the right hand side of 
Eq. (|3.30p has been evaluated in Appendix |B] with the result 



Mji(3ay - 1) — + 



T T \ T — T- 



(3.33) 
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In the case of a three-dimensional system {d = 3), Eq. (|3.33l) reduces to the one previously derived for hard spheres 
In addition, for identical particles previous results tl^il5| obtained for monodisperse gases are also recovered. 
With the right side of Eqs. (|3.30p now determined the algebraic equations 



1) 



T, T, 



rrii + nij 



(3.34) 



can be solved to determine the partial contributions ij^. Their sum then gives the kinetic contribution to the shear 
viscosity, r]'^ . Finally, adding this to the coUisional transfer contribution of (|3.25l ) gives the total shear viscosity. 

C. Heat Flux 

The constitutive equation for the heat flux q*^^^ has contributions proportional to gradients of the densities and the 
temperature, and terms proportional to an applied force (taken to have the same order of magnitude as the gradients) 



q^'^ =-Y.J2 {T^I^i,^3^ Inn, + LyF,) - TAV InT, 

i=i j=i 



(3.35) 



where A is the thermal conductivity coefficient and -Dq,ij are the Dufour coefficients. As in the case of the shear 
viscosity, the transport coefficients -Dq,ij , Lij and A have kinetic and coUisional contributions 



Lij — L^j -\- Lj^j^ A — A " -{- A 



(3.36) 



The coUisional transfer contributions A^, -D^ ^-, and L'ij are given in the preceding paper by Eqs. (7. 14)^(7.16) 



i=i j=i 



^ ^ - (1 + a^j) mj^i^jafjxf^ <{ 2B4 (1 - a^,) (^iy - n 

(2^y - ^iji) pDj 



nrij mimjl ■' 



T, 



rrii ^ mimjl 



lr<T 



(3.37) 



p=l 



2 „j, , , Ti m^fii 



P3 



8B2 
J 

d + 2 



q,pj 



pmi 

T Tl Tfl 

{d + 2) (2p,p-pp,)-^^^Dp 



-2r<T 

ipj 



(3.38) 



s 1 



p=i 



4, + (d + 2)- 



mimp 



-Di 



8B2 

7 

d + 2 



2M: 



■pt Tk 



- {d + 2) {2p,p- fipi) 



mimp 



-D 



P3 



The constants, 5^, are defined by Eq. (jA19p of Appendix [Xl and 



cl 



2B3 
d 



J dvi J dv2/f )(Vi)/f (V2) {gGl + g-' (g • G,,)2 + (1 + ^l,,)g (g • G,,) 



+^J■3^^J■^Jg^ ^^C^^ - [g (g • G,,) +g^]\ , 



(3.39) 



(3.40) 
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Cl„ = ^1^^^! [- (1 - "'f ) (/^^p - 9' + 45 (g • G^p)] ^(Vi)n,a„^/(°)(V2), (3.41) 
where Gip = ^ipVi + /ipiV2. In the first Sonine approximation these integrals have the explicit forms 



X {2/32 + e,e, + (0, + 0,) [{6, + 0,)M^jMr* + + Mj^)]} 



3^(<i-l)/2 



2dT\ 



3/2 



4 (d-l)/2 



(9 In 



7p 



91nn,- 



Oj + Or 



3/2 



- 3 Oj d\TL"fp 
2di + dpdlTinj 



(3.42) 



(3.43) 



where Pip = jiipOp — ^pi9i. With these expressions, the collisional contributions to the heat flux are explicitly known. 
For mechanically equivalent particles, all the above expressions for the collision transfer contributions reduce again 
to the previous results obtained for monocomponent gases [13, [iBl • 



1. Kinetic contributions 
The kinetic parts of the transport coefficients A'^, -D^ i, , if, are defined respectively as 



(3.44) 



(3.45) 



(3.46) 



where Ai(V), S:^(V), and £^^(V) are again the solutions to p.5P " (|3.7p . The leading Sonine approximations of 1 
p.lOp are not adequate to determine the general leading order for these transport coefficients, since the coefficients 
DJ ^ Dij, and Df^ vanish for a simple one component fluid. This would imply a vanishing thermal conductivity A as 
well. Consequently, it is necessary to include here the next (second order) Sonine polynomial 



A(V) ^ /,,m(V) 



(3.47) 



BKV)^ Am(V) 



pn,T, 



: + 2 n,i;3%y^nvj 



(3.48) 



^i(V)^/.,M(V) 



n,T,, 



-YD, 



d + 2 n,Tf 



,S,(V) 



(3.49) 
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where the next order polynomial Si(V) is 

S.(V)= Qm.V^2-^^7^,)v. (3.50) 

In the above equations, it is understood that the transport coefficients Df, Dij, and Dfj are given by Eqs. p.l9p -( 
I3.2ip . respectively. The coefficients A^, , and £ij are the projections along S^: 



A,: 



J dvS,(V). A(V), (3.51) 
J rfvS,(V).Bl(V), (3.52) 



-i J dvS,(V).£|(V). (3.53) 
In terms of Ai, dq^ij and £ij, the transport coefficients A*^, Dq,ij, and L^J become 

i=l 

+ (3.55) 

+ (3.56) 

bmce Df, Dij, and Df, are known from the analysis above it remains to determine Ai, dq^ij and iij. The alge- 
braic equations determining them are obtained by substituting p.47p - (|3.49p into the integral equations p.5P " (|3.7p . 
multiplying Eq. (|3.5p by Si(V) and integrating over the velocity. The results are 

s 

J2il^,-2C^°^S.,)\,^K, (3.57) 



with the inhomogeneities given by 



I] ( 7^£ - lc^°^Sa ) dq,ij = dq,^j (3.58) 



/c=l ^ ^ 



A- = + 2 p n,7;3 ^ _ ^(0)^^^. d{d + 2) n,7;^ 

* 2 Tm^ ^ n,T, ^ 2d WiT 



h^J '^^^^./3(V)/C.,-^[Vv • (V/f )], (3.60) 



2Td 

- (i + 2ninj7;3 /m, ^ cju ~ C^^^^Sa ^ 9C^°\ 1 ainT, 
_ Z^"^^ + ^ — ^ 



(0). 



(3.61) 
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E 



k=l 



nkTk 



D 



ky 



(3.62) 



Use has been made of the results 



y"dvS, •A,(V) 
j dvS,.B^"(V) = - 



d{d + 2) n^Tf 1 



■ ^ / dv5,,^(V)/C.,, ^[Vv • (V/f )], (3.63) 



(i(d + 2) niUjTi dTi 
2 mi 9n,- 



J dvS,(V)-E^ (V) = 



f(o) 



The colHsion frequencies introduced here are 

2 



Hi + 2) „.|l (^E / -ivS. ^ J™I/.,mS,. 4"'] + /dv S. ^ 4°'l/f' ,/..mS.1 



7»j 



d(c? + 2) n,T, 



k\j:jd.s.- j^[h,M^jf] +jd.s.- Ji^'^i/f ,/,mv] 



d(rf + 2) n,T/^ 



(3.64) 
(3.65) 

(3.66) 
(3.67) 
(3.68) 
(3.69) 



The expUcit expressions for the above coUision frequencies are also displayed in Appendix \K\ 

The inhomogeneous terms p.60p and p.6ip involve integrals over the operator ICij.-y. These are evaluated in 
Appendix [B] 



Jdv 5,,^(V)/Cy,^[Vv(V/f)) 



rd/2 



T, 



+ {2d - 5 - 9a^j)^l^J^lJ, + {d-l + 3ay + 60^)/^^^,] + 6— ^^^^^,(1 



(3.70) 



2 Q„ ^,,2 1 



X <^ [(d + 8)^ifj + {7 + 2d- 9a,j)^i^jfi^i + {2 + d + 3afj - 3ay )^|i 
+3fi%{l + oi.jf^ + [(d + 2)^^^- + (2d - 5 - 9a^j)^l^j^lji + {d - 1 + 3a^j + 6al)fi^^^] 



(-- 




-1 









(0)^1 51n7£^ 
9 Inn, 



(3.71) 
(3.72) 



This completely determines the parameters of the integral equations (|3.57p - (|3.59p . Their solution determines A,, 
dq^ij and ^y, and hence the transport coefRcients A*^, D'^ij, and L*"^ through Eqs. (|3.54p - (l3.56p . In the special case 
of mechanically equivalent particles, this description reduces to those previously obtained for a monodisperse gas 

Hlli. 



12 



IV. SUMMARY AND DISCUSSION 

First, a brief summary of the results obtained in the two papers presented here is given. The revised Enskog kinetic 
equation was used as the basis for deriving exact balance equations for the hydrodynamic fields. The pressure, cooling 
rate, mass flux, momentum flux, and energy flux were determined from the kinetic theory as linear combinations of 
the gradients of these fields and an applied force. These are the constitutive equations that convert the exact balance 
equations into closed equations for the hydrodynamic fields (Navier-Stokes equations, at this order of gradients). The 
coefficients in these equations are the pressure, leading order cooling rate, and the transport coefficients. Formally 
exact expressions for these coefficients were obtained in the first paper, expressed as integrals over solutions to integral 
equations via the Chapman-Enskog method for constructing a solution to the kinetic equation. In the present paper, 
these exact expressions were evaluated using approximate solutions to the integral equations. The approximation 
consists of expanding the unknown functions in a complete set of polynomials, and truncating that expansion to 
convert the integral equations to algebraic equations that can be solved by standard matrix methods. Here, the 
detailed forms for those equations or their solutions have been given in terms of the hydrodynamic fields and other 
parameters of the problem. The result is a complete description of the hydrodynamic equations to Navier-Stokes order 
with all parameters determined explicitly from the theory. Due to the complexity of the analysis it may be useful to 
give the specific location of the final results. The pressure is given by Eq. (|2.2p and the lowest order cooling rate by 
Eq. (|2.4p . The equivalence of all species cooling rates determine the species temperatures as well (to this order in the 
gradients). The only transport coefficient for the cooling rate to first order is given by (|2.20|) and the solution to the 
algebraic equations (j2.16|) . The transport coefficients characterizing the mass flux are the solutions to the algebraic 
equations p.l9p - p.2ip . The momentum flux has two transport coefficients, the shear and bulk viscosities. These 
are given by Eqs. (|3.23p - (|3.26|) and the solutions to the algebraic equations (|3.34p . Finally, the transport coefficients 
for the heat flux are given by Eqs. (|3.36p - p.39p . (|3.54P " (|3.56p . and the solutions to the algebraic equations p.57p - 
p.59p . Detailed forms for the collision frequencies and other input functions are presented in the Appendices. These 
explicit "constitutive relations" together with the exact macroscopic balance equations for species mass densities, flow 
velocity, and temperature (Eqs. (4.12), (4.14), and (4.24) of the previous paper) complete the practical description of 
Navier-Stokes hydrodynamics for a moderately dense, multi-component granular mixture based on the revised Enskog 
kinetic equation. 

In the remainder of this Discussion the relationship to several previous works on kinetic-theory-based descriptions of 
granular mixtures in the literature is considered. These contributions, which are abbreviated according to the initials of 
authors and the flnal two digits of the publication year, are listed in Tables HI and [III To help put the current effort in the 
context of previous works. Tables HI and [Til also list the applicability of each contribution (dimensionality and number of 
species), the starting kinetic equation (Boltzmann or Enskog), and the specific mechanics and assumptions used in the 
derivation process: the standard Enskog theory (SET) vs. the revised Enskog theory (RET), pair correlation function 

Xij\ solution method, order and base state of expansion, single particle velocity distribution function [Maxwellian (M) 
and non-Maxwellian (nM)], energy distribution [energy equipartition (EE) and non-energy equipartition (nEE)], and 
hydrodynamic variables. In the subsections below, each of these various treatments are detailed and their implications 
are discussed. 

Dimensionality 

The number contained in Table U refers to the dimensionality of the particles. Namely, 2D refers to circular 
particles (disks) that are constrained to motion in a plane, while 3D refers to spherical particles that can move 
in all three dimensions. Although systems of practical importance are three-dimensional, two-dimensional theories 
have appeared in the literature for purposes of comparing with molecular-dynamics (MD) simulations. Specifically, 
early MD simulations were performed in two dimensions due to computational constraints, whereas three-dimensional 
simulations are now common. The WA99 theory is applicable to 2D only, the JM87 and GHD theories are derived 
for both 2D and 3D, while all remaining theories JM89, Z95, H06, GD02, R03, IA05, and S06 are for 3D systems. 

Number of Species 

Most mixture theories to date have been targeted at binary mixtures, in which the two species can differ in size, 
mass, restitution coefficient, and density (JM87, JM89, WA99, HOI, GD02, R03, and S06). Some recent theories have 
been derived for a more general system of s distinct species (Z95, IA05, and GHD). 

Kinetic Theory 

Two kinetic theories have been used in the development of hydrodynamic equations for mixtures, namely the 
Boltzmann equation and the Enskog equation for hard spheres. The difference between these two equations stems 
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from the treatment of the two-particle distribution function fij. For the Boltzmann equation, is assumed equal to 
the product of the two single-particle distribution functions (/^ — fifj)- This lack of spatial and pre-coUisional velocity 
correlations between the two particles restricts the Boltzmann equation to dilute systems. The Enskog equation, on 
the other hand, accounts for positional correlations (but not pre-collisional velocity correlations) via the equilibrium 
pair correlation function at contact Xy, namely fij = Xijfifj (the Enskog approximation). More specifically, Xij 
accounts for excluded volume effects encountered in denser flows, and thus the corresponding Enskog kinetic theory 
is applicable to moderately dense flows. The Enskog approximation is expected to deteriorate at higher densities as 
ring collisions and their associated velocity correlations become important. The prevalence of such correlations has 
also been found to depend on the restitution coefficient; namely, velocity correlations have been observed to increase 
as the restitution coefficient decreases. Thus, theories based on the Enskog equation are applicable to both dilute and 
moderately dense flows (JM87, JM89, Z95, WA99, HOI, R03, 1A05, and GHD), while theories using the Boltzmann 
equation as a starting point are restricted to dilute flows (GD02 and S06). 

From a practical perspective, the upper limit of concentration that a given kinetic theory should be applied to 
depends on the desired level of accuracy. As a quick gauge to the range of validity of the Boltzmann equation, the 
limiting case of a monodisperse system is considered using the theories of Garzo and Dufty [Ti] and Lutsko [3], 
which are based on the Enskog equation and thus are applicable to both dilute and dense flows. According to these 
theories, the dense coUisional contributions to the pressure (absent in the Boltzmann equation) are 4% and 18% for 
solids volume fractions of 0.01 and 0.05, respectively, at a = 0.9. A similar estimate of the range of validity of the 
Enskog equation is not available, since the impact of velocity correlations on granular hydrodynamics has not been 
extensively studied. 

SET vs. RET 

For those contributions listed in Table U that employ the Enskog equation, two different approaches are possible - 
the SET and the RET. As mentioned previously, the difference between SET and RET traces to the choice of the 
pair correlation function, Xij^ ■ SET, xlj'' is a function of concentration (i.e., depends on local value only) at a 
single position of interest, whereas for RET Xij^ is treated as a functional of concentration (i.e., depends on the local 

value and its gradient) at the two particle centers. In SET, the location (e.g., midpoint) at which to evaluate xfj^ for 
mixtures is unclear [l^. Regardless of the choice, however, the resulting diffusion force is found to be inconsistent 
with irreversible thermodynamics, unlike in RET [l9|. 

The implications of the SET vs. RET treatment on the resulting theory depends on the type system being 
examined. For the case of monodisperse systems, SET and RET lead to the same Navier-Stokes transport coefficients 
0, l20|, \M ll but different inhomo geneous equilibrium states [l^, |20|, |2l| . For mixtures, however, different Navier- 
Stokes transport coefficients are obtained from SET and RET. More specifically, although the fluxes appearing in 
the momentum and energy balances are the same for SET and RET [q and P-y^ in Eqs. (4.13) and (4.14) in the 
companion paper], the diffusion flux [joi in Eq. (4.12) in companion paper] takes on different forms [l^. Although the 
quantitative impact of such differences on segregation predictions has not been investigated in detail, it is clear that 
RET is the appropriate approach since it is consistent with irreversible thermodynamics. The Enskog-based theories 
of JM89, Z95, WA99, and GHD use RET, whereas JM87, HOI, R03, and IA05 utihze SET. 

Pair Correlation Function 

As noted in TablelH various forms of the pair correlation function Xtj"^ have been used in conjunction with mixtures. 
Its explicit forms can be found in the references displayed in the table. 

Solution Method 

With the exception of Zamankhan [ij] , all of the efforts to date on mixtures have implemented a Chapman-Enskog 
(CE) expansion to solve the kinetic equation. By deflnition, the CE expansion involves a perturbative expansion 
about low Knudsen numbers (where Kn is defined as the ratio of the mean free path to the characteristic length of 
the mean-flow gradients) or "small gradients," and thus is not applicable to systems in which free-molecular (non- 
continuum) effects play a non-negligible role. The Grad moment method, as was employed by Zamankhan [4j] does 
not contain similar restrictions, though the derivation is necessarily more complex and thus has not been performed 
without resorting to other simplifying assumptions (e.g., equipartition of energy). 



Order of Chapman-Enskog Expansion 
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For each of the entries in Table |T] that employ the CE expansion to solve the kinetic theory equation, all expansions 
are carried out to first order in spatial gradients - i.e., to the Navier Stokes (NS) order. Nonetheless, evidence of 
higher-order effects has been noted in a range of granular flows, and Burnett-order effects (second order in spatial 
gradients) in particular have been shown to be linked to the anisotropy in the stress tensor '24'|. Examples of systems 
in which significant higher-order effects have been identified include the upper region of an open-ended, vibro-fluidized 
bed [2^ [2^ , the vanishing heap of a vibr ated, granular material dilute flow around an immersed cylinder d, 
simple shear flow [295 , the Knudsen layer a djace nt to a thermal (energy-providing) boundary [30| , and in the continuum 
interior of a thermally-driven granular gas [3l|. For the case of monodisperse systems, several theories[24, 32, 33, 3j| 
and boundary conditions [2^, [2a [s^ have been developed that account for such higher order effects in a limited class 
of systems (dilute, sheared flow, open boundary, etc.). Analogous work has not been reported for mixtures [SGllSTj. 
however, and thus the incorporation of higher-order effects does not serve as a differentiator between those theories 
listed in Table HIl that employ the CE expansion. Nonetheless, the possible presence of such higher-order effects should 
be kept in mind when comparing these theories with experiments and simulations for purposes of validation. 

Base State of Chapman- Enskog Expansion 

As described in the companion paper, hydrodynamics results from a "normal solution" to the kinetic equation, 
whose space and time dependence occurs only through the hydrodynamic fields. The CE expansion is a systematic 
method for constructing this normal solution as an expansion in powers of the Knudsen number, or spatial gradients 
of the fields. At zeroth order in these fields, the kinetic theory determines the form of the distribution function 
to be the "local equilibrium" Maxwellian for molecular fiuids. However, in the presence of dissipation the kinetic 
theory requires a different solution at zeroth order, the "local homogeneous cooling" (HCS) distribution. The HCS 
distribution agrees with the local Maxwellian only for a = 1. It is possible to make an expansion in both Knudsen 
number and (1 — a), in which case the lowest order term is indeed the Maxwellian. Such a double expansion is 
necessarily limited to asymptotically weak dissipation. The theories of JM87, JM89, WA99, HOI, R03, IA05, and 
S06 are of this type with the corresponding implicit limitation. The theories of GD02 and GHD, on the other hand, 
expand only in the Knudsen number, with HCS as the leading order solution and hence no a priori limitation on 
the degree of dissipation. Quantitatively, the difference between the two types of expansions has been examined in 
monodisperse systems via a comparison of the dissipation rate obtained in MD simulations [Slj. As an example, at 
a volume fraction of 0.3 and a = 0.75, an error of 23% is obtained when using a theory [sst based on an expansion 
about a — 1 , while an error of 7% is obtained when using a theory [l3| based on an expansion about the HCS. For 
both theories, the level of mismatch is found to increase with concentration and dissipation levels. 

Order of Sonine polynomial expansion 

To make the analytical evaluation of the collision integrals possible, a truncated Sonine polynomial expansion is 
employed. All theories listed in Table [J employ the lowest, non-zero order of the polynomial expansion (leading 
term), except for the weak dissipation theory of S06 which carries out the expansion to third order. Nonetheless, the 
transport coefficients of S06 agree with those of GD02 in the common domain of validity, namely in the nearly elastic 
limit 1^. 

The accuracy of a given Sonine polynomial approximation can be tested via comparison with discrete simulation 
Monte Carlo (DSMC) results. DSMC provides a numerical solution of the starting kinetic (Boltzmann or Enskog) 
equation for a specific system, and thus provides a check for the existence of a normal solution, the order of the 
gradient expansion (e.g., Navier Stokes), and the truncated Sonine polynomial. For dilute granular mixtures, good 
agreement between DSMC and the theory of Garzo and Dufty [ll| was found for the shear viscosity, even for strong 
dissipation [3^. This agreement for the shear viscosity is also present for moderately dense systems [l^- A similar 
comparison was carried out for the diffusion coefficient for a system with an impurity (40| . The agreement is again 
excellent, except when the size or mass disparity is large, in which case the second Sonine approximation leads to a 
significant improvement over the first Sonine approximation. 

Single Particle Velocity Distribution 

In an effort to simplify the evaluation of collision integrals, some previous works have assumed that the single particle 
velocity distribution function is Maxwellian (M) . Strictly speaking, such an assumption is valid only for perfectly elastic 
spheres in equilibrium [4l| . Numerous experimental, theoretical, and simulation studies of inelastic, monodisperse 
systems have indicated that the distribution function departs from Maxwellian [H, [H, [H, [H, [H, \^ [H, \^ . For the 
case of perfectly elastic mixtures, an estimate of the impact of this effect is given by Willits and Arnarson Q, who 
compare the shear viscosity predictions of two theories with that of MD simulations. For this particular system, the 
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only effective difference between the two theories is that one contains a Maxwellian assumption and the other does 
not. For disks with a diameter ratio of 1.25 and over a range of solids fractions from 0.05 to 0.4, the non-Maxwellian 
(nM) theory exhibits very good agreement with the simulations. The Maxwellian-based predictions, on the other 
hand, are significantly lower in value since the Maxwellian assumption precludes the kinetic contribution to the shear 
stress and simplifies its collisional counterpart. Although similar studies on the impact of non-Maxwellian effects on 
other constitutive quantities are not available, it is clear that a non-Maxwellian treatment is critical for the accurate 
prediction of shear stress. Of the mixture theories presented in Table HI those of JM87, HOI, and R03 are based on a 
Maxwellian assumption, while those of JM89, WA99, Z95, GD02, S06, and GHD incorporate non-Maxwellian effects. 
The theory of IA05 takes a hybrid approach, in which collision integrals involving unlike particles are evaluated based 
on Maxwellian distributions, whereas collisions integrals involving like particles account for non-Maxwellian effects. 



Energy Distribution 



Similar to the assumption of a Maxwellian velocity distribution, the assumption of an equipartition of energy (EE) 
has also been made periodically in an effort to simplify the evaluation of collision integrals. Again, an equipartition 
of e nerg y between unlike particles is only ex pec ted for a perfectly elastic system in equilibrium. Numerous theoretical 
[H, l39l. Isoj , simulation [5ll, [H, [13, [H, [H, , and experimental [13, [IE [111 studies of granular materials provide 
evidence that an equipartition of energy does not exist, and that the level of non-equipartition of energy (nEE) increases 
as the restitution coefficient decreases (dissipation increases) and as the mass ratio gets further from unity. The impact 
of such non-equipartition has recently been evaluated in the context of species segregation. In particular, Galvin, 
Dahl and Hrenya [63| have found that the driving forces for species segregation that arise from non-equipartition are 
significant over a moderate range of parameters for the case of a thermally-driven system. Furthermore, for the case 
of an intruder particle in the presence of gravity, both Brey et al. 61 1, Garzo ^62*1, and Yoon and Jenkins 'gs'I have 



found that the direction of species segregation may reverse due to non-equipartition effects. Together, these studies 
indicate the importance of non-equipartition in a variety of systems. The mixture theories of JM87, HOI, GD02, R03, 
IA05, and GHD include the effects of non-equipartition, while those of JM89, WA99, Z95, and S06 do not account 
for non-equipartition. 



Hydrodynamic Variables 

The appropriate choice of hydrodynamic variables to include in a hydrodynamic theory depends on the timescale 
associated with a given variable [641] . First consider the case of a molecular gas, in which two timescales are relevant 
- the "kinetic" time scale and the "hydrodynamic" time scale. The kinetic time scale is "fast" and representative of 
the time between collisions (i.e., mean free time). The hydrodynamic time scale is "slow" and is set by the gradients 
in the system; this is the timescale over which macroscopic variables change. 

It is important to note that the velocity distribution function fi (and thus the Boltzmann and Enskog equations) 
carries information with it on the kinetic time scale. However, the macroscopic variables (n^, U, and T for the case of 
a molecular gas) correspond to velocity moments of fi that are collisional invariants - i.e., Ui, U, and T are conserved 
upon collision. As a result, they will stay constant over a time on the order of a mean free path, and vary in time 
over the much longer hydrodynamic time scale. It is this property which defines the appropriate (minimum set of) 
hydrodynamic variables needed to describe a system. 

To illustrate the above point, consider the time evolution of mixture temperature (T) and species temperature (T^), 
again for the case of a molecular gas with more than one species present: 

DT 

— = (energy flux) 

DTi ^ a \ ( collisional exchange (source/sink) \ 

Dt ~ ^^'^^^^•^ ^ between unlike species j 

Note that T is a collisional invariant (the total energy is conserved upon collision and thus no "source" term is 
present), but that Ti is not a collisional invariant since the energy associated with a given species can change upon 
collision with another species. Also note that the time scale of the flux term is set by gradients (e.g., Fourier's law), 
whereas the time scale of the source term is collisional in nature and thus relatively "fast." Hence, the mixture 
temperature T is characterized by one (slow) timescale and the species temperature is characterized by two (slow 
and fast) timescales. The corresponding physical picture is the following: each species is expected to have a rapid 
relaxation after a few collisions to a distribution near local equilibrium (local mixture temperature), and then all 
temperatures evolve according to the slow hydrodynamic mode. Thus, the relevant hydrodynamic (macroscopic) 
variable is T since it is associated with a conserved quantity; all other variables (T^) are enslaved to it. Thus, solving 
a separate balance for Ti is superfluous - indeed this is not done for molecular gases. 
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Next, consider a granular system in which the particles engage in dissipative collisions. Three timescales need to 
be considered, as illustrated by the balances below: 

(inelastic \ 
dissipation 
(sink) J 

, n \ i \ / collisional exchange (source/sink) \ 

(energy Hux) + dissipation + i ^ i-i ■ 

^ 'I (sink) / ^ between unlike species ) 

As with the molecular systems, the energy flux is characterized by a slow (hydrodynamic) time scale and the collisional 
exchange between species is characterized by a fast (kinetic) time scale. In the granular system, inelastic dissipation 
also occurs, and its corresponding time scale is not as obvious. More specifically, the dissipation rate depends not 
only on the collisional frequency (kinetic time scale), but also on the value of restitution coefficient. As a approaches 
unity, the dissipation rate goes to zero, and the time scale becomes large (hydrodynamic time scale). In other words, 
it is unclear a 'priori whether this term is characterized by a fast or slow timescale. 

The previous observation leads to the question: for moderate (smaller) values of a, does the time scale become fast 
enough to approach the kinetic time scale or is there still a separation of time scales? If the latter is assumed (that 
the time scale associated with inelastic dissipation is much longer than that of the kinetic time scale), Garzo and 
Dufty [l3] have shown that in a HCS two-component system, (i) the cooling rates associated with mixture and species 
temperatures are equal and (ii) the ratio of species temperature (T1/T2) remains constant - i.e., all temperatures 
decay at the same rate but their ratios remain constant. Further note that T1/T2 is not equal to one for unlike 
particles - i.e., a non-equipartition of energy is predicted. To test this assumption of the separation of timescales, 
MD simulations of a binary mixture in HCS were performed by Dahl et al. The simulation results indicate that 
a constant value of T1/T2 is achieved after only a few collisions, and this behavior was confirmed over a wide range 
of parameters. 

As a result of the aforementioned finding that the timescale associated with inelastic dissipation is hydrodynamic 
(slow) in nature, the physical picture is analogous to that of a molecular system Namely, each species has a rapid 
relaxation to a local "equilibrium" state (now characterized by a constant value of non-equipartition) - HCS - and 
is then enslaved to slow (hydrodynamic) evolution. Thus, the relevant hydrodynamic variables are the same - n^, U, 
and T - and the balance of additional variables (such as Ti) would be superfluous. 

The practical implications of using ni, U, and T as the hydrodynamic variables in granular flows are twofold. 
First, only a single balance for T is needed, instead of separate balance for each Ti. The reduction in the number of 
governing equations is expected to lead to a considerable decrease in computational overhead. Second, the level of 
non-equipartition is in a state of local "equilibrium" due to the fast time time scale of Ti, and thus T1/T2 depends 
on local values of flow field variables and particle properties (mixture composition, a, etc.). Hence, even though 
a non-equipartition of energy is indeed present, it does not appear explicitly in the transport coefficients. Instead, 
its dependency on the flow field variables is incorporated into the transport coefficients - i.e., the effect of non- 
equipartition is implicitly contained in the transport coefficients. It can be solved for explicitly once the flow field 
variables have been solved for using, for example, the relation derived by Garzo and Dufty (l^ . 

Of the non-equipartition theories listed in Table [Hi GD02 and GHD use T as the hydrodynamic variable for the 
granular energy field, whereas the other non-equipartition theories use Ti (JM87, HOI, R03, IA05). The choice of 
hydrodynamic variable(s) is a non-issue for theories which invoke an equipartition assumption (JM89, Z95, WA99, 
S06) since they inherently assume T = Ti and thus use T as the energy variable. (It is worthwhile to note that S06 
defines temperature differently than others, namely Ts — 3T. Furthermore, IA06 uses a species temperature TijA, 
which is defined in terms of velocity fluctuations relative the mean species velocity, rather than the mass-averaged 
velocity used for T.) 

Recapitulation 

Unlike previous Enskog-based (dense) theories, the current effort is based on an expansion about the homogeneous 
cooling state and employs rij, U, and T as the hydrodynamic variables. The former extends the range of validity to 
strong dissipation levels while the latter results in fewei^balance equations, thereby reducing the computational cost. 
This approach was first put forth by Garzo and Dufty 11] who instead used the Boltzmann equation as their starting 
kinetic equation; the current work thus extends the domain of applicability from dilute to moderately dense flows. 



DT 
iDt 

DT, 
dt 
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TABLE I: Hydrodynamic Descriptions of Granular Mixtures 
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(0) 
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JM87 


2D and 3D 


2 


Enskog 


SET 


2D f2l 3D [651 


[3] and [5] 


JM89 


3D 


2 


Enskog 


RET 


3D [65] 


[4] 


Z95 


3D 
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Enskog 


RET 


[66] 


[6] and \67} 


WA99 


2D 


2 


Enskog 


RET 


2D [2] 


m 


HOI 


3D 


2 


Enskog 


RET 


[65] 




GD02 


3D 


2 


Boltzmann 






m 


R03 


3D 


2 


Enskog 


SET 


[68] 


[9] 


IA05 


3D 


s 


Enskog 


SET 


m 


[10] 


S06 


3D 


2 


Boltzmann 






current work 


GHD 


2D and 3D 


s 


Enskog 


RET 


2D \70}, 3D [71J 



TABLE II: Hydrodynamic Descriptions of Granular Mixtures 
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Hydro- 
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Expansion 
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Expansion 


Velocity 


Energy 
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Ref. 


Method 


Order 


Base State 


Order 
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Variables 


JM87 


CE 


NS 


a = 1 


1st 


M 


nEE 


rii, Ui, Ti 


JM89 


CE 


NS 


Q = 1 


1st 


nM 


EE 


U„ T 


Z95 


Grad's method 








nM 


EE 


13 moments 


WA99 


CE 


NS 


a = 1 


1st 


nM 


EE 


n„ U,, T 


HOI 


CE 


NS 


a = 1 


1st 


M 


nEE 


Hi, Ui, Ti 


GD02 


CE 


NS 


HCS 


1st 


nM 


nEE 


n^, U,, T 


R03 


CE 


NS 


Q = 1 


1st 


M 


nEE 


Hi, Ui, Ti 


IA05 


CE 


NS 


a = 1 


1st 


M/nM 


nEE 


Hi, JJi, TijA 


S06 


CE 


NS 


Q = 1 


3rd 


nM 


EE 


Hi, Ui, Ts 


GHD 


CE 


NS 


HCS 


1st 


nM 


nEE 


n„ U,, T 
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APPENDIX A: COLLISION FREQUENCIES 

In this Appendix we display the expressions for the coUision frequencies appearing in the evaluation of the kinetic 
contributions to the transport coefficients and in the first order contributions to the cooling rate. As noted in the 
main text, most of these frequencies (those corresponding to the transport coefhcients) have been obtained in the low- 
density limit (except for the factors ^ ) [3 for arbitrary dimensions by considering the Maxwellian approximation 
for the zeroth-order distribution functions fj^^^. Their expressions are 



^ a-i (0) _^ sfO^ + e, 

V 2 / j^i ^ * 



1/2 



(Al) 



dr(i) 



OiOj 



1/2 
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r(f) d(d + 2)' 
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In the above equations, vq = y^2T/m is the thermal velocity, 6i = miT/mTi, and we have introduced the dimensionless 
quantities 

A,, = id + 2)(2/3y + 6,) + ti,,{e, + 9,) {{d + 2){l~ a,,) - [(11 + d)a,, -hd- 7]A,0ri} 

+3(d + 3)4.0ri + 2^l% {2al - +d+l^ Q-^O, + e,f - {d + 2)e,er\e. + 9,), 



(A9) 



Cy ^ {d + 2)(2/3,, - e,) + fi,,i9, + 9,) {{d + 2)(1 - a,,) + [(11 + d)ai2 -5d- 7]/3„-0-i} 



3 



(AlO) 



E,, = 2Ai2^0-2(^^ ^ ^2al - + d + l) [(d + 2)9, + {d + 5)9,] 

-Mj.(0, + 9,) {(3,j9^^[{d + 2)9, + {d + b)9j][{ll + d)a,, -5d-7] 

-9j9r^[20 + d(15 - 7a,j) + d\l - a,j) - 28a,j] - {d + 2)\l - a,,)} 

+3{d + 3)/32 6'r2[(d + 2)9, + {d + 5)9,] + 2p,,9r\{d + 2)26*, + (24 + lid + ^2)617 

+ {d + 2)9j9r\d + 8)61, + {d + 3)07 - (d + 2)(6i, + 9j)9r'^9,[{d + 2)6*,; + (d + S)^^, 



F„ = 2^%9J\9, + (^24 - ^a,, + d + 1^ [(d + 5)0, + (d + 2)0,] 

-/i,,(0, + 0,) {Aj^7'[(rf + 5)0, + (d + 2)0,][(11 + d)a„ - 5d - 7] 

+0,07i[2O + d(15 - 7a,j) + d^{l- a^- ) - 28a,j] + (d + 2)^(1 - a^- )} 

+3(d + 3)/32 6i72[(rf + 5)6)^ + + 2)0j] - 2/3,j9j^[{24 + lid + d2)0, + (d + 2)^0^ 

+ (d + 2)0,071 [(d + 3)0, + (d + 8)0j] - (d + 2)(0, + 9,)9j\{d + 3)0, + (d + 2)07, 



(All) 



(A12) 



where = ^ij0j - /ij,0,. 

It only remains to evaluate the collision frequencies tpij needed to get the first order contribution ^(^-i) to the 
cooling rate. These frequencies have not been previously determined in the Boltzmann limit. To compute them, we 
use the property 

ydvi/i(vi)j(^[Vi|<i>„$7 = xf'yff' J dvi J dV2<i>,(vi)<i>j(v2) 

X y dS e(^ • g)(S • g) [h{V[) - h{Vi)] , (A13) 

with 

V'i=Vi-Aij,(l + a,j)(5^-g)5^. (A14) 
To determine ipij, let us consider first the integral 



I^-l d-(^) 4"^[/rJ.MF,]. (A15) 



Use of Eqs. (|A13p and (|A14p gives 



U^xf4~'^ J / 'iv^/f (Vi)/,,m(V2)f,(V2) J daeia-g){a-g) {v;'-v,') . (ai6) 
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The collision rule (|A14p yields 



+ a,,){a ■ ■ Vi) [V,' + 4(1 + a,,f{a ■ g)^] 

The integration over the solid angle in (|A16[) leads to 

4i?3 2 



d + 3 



+ a»j)ff(Vi • g) [(d + 3)V^^ + 4^(1 + ay) V] 



where 



i3fc = y dse(^ g) •g)'= 



_ ^(rf-l)/2 



Therefore, the integral (jA16p can be written as 



lib — 



_ (0) 



d + 3 



where we have taken the Maxwellian form for fj^^"^ and have introduced the dimensionless integral 



X \ 4^jj(l + a,j) 



3(ci.x)^ + ^5^c? + M,^.(l + «.,)V 



-4(ci • x) [[d + i)cl + An%{l + a,jfx^\ } , 



and the operator 



d d{d + 2) 



d0i 



(A17) 



(A18) 



(A19) 



(A20) 



(A21) 



(A22) 



In addition, = Vi/wq and x = g/fo- The integral (|A2ip can be performed by the change of variables {ci,C2} 
{x,y} where y = 9iCi + 0jC2 and the Jacobian is {Oi + 0j)~'^- With this change the integrations can be done quite 
efficiently by using a computer package of symbolic calculations. A lengthy and careful algebra gives 



X {-29^ [{d + 3)9, + {d + 2)9,] + fi,, (1 + a,j)9, {9, + 9,) 
d"^ + 5d + 6, ' 



(11 + d)9, 



d + 3 



-84 (1 + a.,,)' 9, [9, + 9,f + 2^i% (1 + a,,)' {O; + Ojf} 
Use of (|X20)) in ((J23l) leads to the final expression for : 



(A23) 



l%h — 



(0) d-\ /I , Nfl3/2^-3/2 



4r(C 

X {[{d - l)9j + {d + 2)9,] [29, + 3fi,, (1 + ay) (9, + 9j)] 
-244 (1 + ay)' {9, + 9j)' + (1 + ayf 9j^ {9, + 9,)'} 



(A24) 
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Following similar mathematical steps as made before for one obtains 

X [20, + ( 1 + ay ) {9, +e,)]. (A25) 
Combining Eqs. (jA24p and (|A25p . one gets the result 

X {[{d - 1)9, + {d + 2)9^] [26, + 3mj, (1 + ay) {9, + 9,)] 

-24fi% (1 + a^jf (9, + 9,f + 30m,', (1 + OJ^ + 0,)' 

-{d+2) {9, + 9j) [29 J + (1 + ay) (0, + 0,)]} . (A26) 

The remaining integrals needed to determine the collision frequencies ij^a and V'y" can be also obtained by performing 
identical mathematical steps. Their expressions are 

dvF,{y)jf [U,mF,, /]"'] = !^^„,n,x!?<riz;„^^, (1 + ay) {9,9,)-^/\9, + 9,)-^'^ 

X {-2 [(45 + 15^)6*1 + 3(38 + 13d)6»i6'f 
+8(11 + Ad)9}9j + 8(2 + d)9l] 

VdyLj^ (1 + ay) (6*, + 9j) [(55 + M)9] + 9(10 + d)6',6i, + 4(8 + d)9^] 
-24^i% (1 + ay)' {9^ + 9,f (59, + 49,) + 30^i% (1 + ay)' (9, + 9,)^ 
+ {d + 2)9, {9, + 9j) [2(40, + 30,) - 3^*,, (1 + ay) [9, + 0,)]} , 



(A27) 



d-vt^{\)J^^ ,/j,j\/-f'iJ - ^^d^ "-tX,, (1 + aij)^'o0,; 



X [I0a| + 22a2 + 11a,, - 3] , (A28) 

j d.F,iy)J^[kMFUr] = ~^<^-ri^ + o..^v,9-''' 

X [SOafi - 126a-, + lllau + 16(i(3a„ - 7) - 137] . 

(A29) 

The expressions for the frequencies ipu and ■i/'y can be easily obtained from Eqs. (|A27p - (|A30p when one takes into 
account their definitions (|2.17p and (|2.18[) . 

In the case of mechanically equivalent particles (x,-"'' = x'"'^\ '^i = ^^nd ay — a), Eqs. (IA28|) and (|A29p yield 



dv^^(V)(j(°)[/(°),/Mf^(V)] + j(°)[/MF(V),/(")]) = !L__^(o)^.-i„^,^.. ^ 

X (30a' - 30a2 + 105a + 24da - 56d - 73) 



(A30) 

This expression coincides with the one previously derived for a monodisperse granular gas [l5l |. 



APPENDIX B: COLLISION INTEGRALS 



In this Appendix, we provide some of the mathematical steps to compute the different collision integrals involving 
the operator /Cy ,y[A]. This operator is defined as 



/Cy,^[A,] = crf^xf j dv2 J d^e(^-gi2)(?-gi2)^?. 



«,^Vr(V'/)X,(V^') + /f (Vi)A,(V2) 



(Bl) 
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To simplify all these type of integrals, we use the property 

J dvi/i(Vi)/C.,-^[X,(V2)] = -xf4j dvi J dv2/f)(Vi)X,(V2) 

X J daQ{a-g){a-g)d^[h{v[)-h{Vi)] 



(B2) 



where V'^ is defined by Eq. (|A14p . 

Let us start with the collision integrals appearing in the evaluation of the mass flux. One of them is 

Id = J dvim,yi,^/C„,^[Vv. •(V2/f^)]. 

Use of the identity (|B2| in (lB3| gives 

lD=B,af^xfm,ti,,{l + a,,) J dV, J rfV^/f (Vi)Vv. • (Vs/f (V2)) 
The integral ljB4[) can be exactly evaluated and the result is 

Id = -2dB2ninjafjxfj^ + 



(B3) 



(B4) 



(B5) 



The remaining integrals corresponding to the mass flux can be computed by using similar mathematical steps as those 
made before for Id- The results are 



J \ J 



(B6) 





{-- 








mi J 



(B7) 



njTTii dlnrij 

The collision integral involved in the evaluation of the pressure tensor is of the form 

I,^ ^ J dvim^Vi^^Vi^i3lCij^^[dvfJj°^] 

= j rfvi j dv2/.(Vi)(ay,,,/f (V2)) 

X j da Q{a ■ g)(^ • g)?„ (K7K/3 - yi,iyi,p) - 

where the identity (jB2p has been used. The scattering rule (|A14|) gives 

^,7^/3 ~ ^i,7^i,/3 = -Mij(l + ■ g) [Gy,-y(T/3 + Gij^pa^ + ^j.,{g^ap + gpd^ 

+ a,j){a ■ g)d^df3\, (B9) 

where Gy = Vi + fJ,jiV2. Substitution of Eq. (|B9p into jBSj allows the angular integral to be performed with the 
result 



(B8) 



J da e{a • g)(a • g)a^ (K7K/3 - VmVi.p) 



Bo 



d + 2 



fiji{l + a^j) [id + 3)g'^Gi-j,p 



+ fj.ji{l + d- ia^j)g'^gfi + 2(g • Gjj)^^] . 



(BIO) 
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With this result the integral becomes 



dB2 



-Hi 



n,xf4'n,t^i^^l + a,,) I dv, I dv2/r(Vi)/f ^(V2) 



X + 2)^,,(3ay - - 4(d+ 2)(g • Gy)] 



fJ-ji 



Tr T, 



rrij nin 



T- — T 
TOi + rrij 



(Bll) 



From Eq. (|B11|1 it is easy to get the expression (|3.34p . 

To evaluate the collision integrals appearing in the determination of the heat flux one needs the partial results 



s.(v;) - S,(Vi) = ^(1 + a^,)^i,,i^ ■ g) { [(1 - • g)' - Gl - 45 



-2^ijiig ■ Gy ) + 2(1 + a,J)^lJ,{a ■ g){a ■ G.j) + {d + 2) 
- [(1 - ■ g) + 2(ct • Gy)] (Gy + Aij^g)} , 



T, 



(B12) 



d^e(^ •g)(a.g) ^ • [S,(Vi)-S,(Vi 

1 



3 Bo 



;(l + ay)/^j. 



2d + 2' 



-A^f, [(d + 2)-3ay(l 



^/i,,(7 + 2d-9a,,).g2(g.Gy) + ^(g.G,^2 V" - „2 



(d + 2)2 T, 



(B13) 



The corresponding integrals associated with the heat flux can be explicitly evaluated by using Eqs. (|B12[) and (|B13[) 
and the same mathematical steps as before. After a lengthy algebra, one gets the expressions p.70p . p.7ip . and 

Let us consider now the integral appearing in the evaluation of the cooling rate. To do that, we compute first the 
collision integral 



(B14) 



In order to evaluate it we use Eq. (jA17[) and the angular integrations 

3^2 



J d^e(^ •g)(^-g)a^ 



(d + 2)(d + 4f-(^ + "-) 
X + a,,) [2(Vi • g)^g0 + 2(Vi • g)^^!.^ 

1 



-,{d + 6)g'V,'gp 



--{d + 4)V,^2{V,-g)g0+g^Vi^p] 
-44(f + a^jfg' [4(Vi • g)gp + g'Vi,p] } . 



(B15) 
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With these results, the integral becomes 



X {^,,(1 + ay) [2(d + 4).gVi2 + 8(Vi • g)^ + 5^(1 + ay ) V - 16Mj.(1 + a«j)5'(Vi • g)] 



--{d + 2)V,\V,-g) 



^2^ 



T 



3dB2ninjXij + ay) i [//^^(l + ay) (16 + 2(i - 16^-,i(l + ay) 



+5M|,(l + «y)') --(d + 2) 



+ 5/ij',(l + ay)'3 



77 



2 T^Tj 
(i + 2 mim 



-flj^i^ + ay) [rf(<i + 4) + 4 - 8(d + 2)^ji{l + ay) + 5(d + 2)^2^(1 + ay)^] I 
j ) 



(B16) 



The final expression (|2.19p can be easily obtained from Eqs. (jBlip and (|B16p . 

Finally, note that the integrals Cy and C'^j defined by Eqs. (|3.40p and (|3.4ip . respectively, can be easily computed 
by using the change of variables written below Eq. (IA22P . After some algebra and using the Maxwellian approach for 
the distributions one gets the results (|3.42l) and (|3.43|) . 



APPENDIX C: CHOICE OF xif AND luj 

This Appendix deals with the choice of the pair correlation function Xij^ and the functional derivative I^j . A good 
approximation for Xij^ in two dimensions (d = 2) is given by [t^ 

l-^ 16 (1-0)2 ay 16 0(1 -^)\a,j J ' ^ ' 

where (p = J2i=i 4'i is the total solid volume fraction, (pi = niirul /A is the species volume fraction of component i and 
(3 = 7i'(X]i=i "iO'i)/4- In the case of hard-spheres {d = 3) we take the following approximation for x|°^ [Tlj 

1-0^2(1- 0)2 ay ^ 2 (1 - 0)3 ay ) ' ^ ' 

where now 0^ = riiTraf /6 and /? = 7i'(X]i=i ''^i'^I)/^- 

The parameter /^f j is chosen to recover the results derived by Lopez de Haro et al. [73| for ordinary polydisperse 
mixtures in the context of the RET. To do that, for the sake of simplicity, we assume that the temperature, the 
pressure and the flow velocity are homogeneous so that only the spatial gradients associated with the partial densities 
will be considered. In this simple case, for elastic collisions (ay = 1), the first-order distribution function is given by 

f'f'^ = Yl,]=i ' Vlnnj where B^(V) verifies the integral equation 

(CB^),^^^ = -V^n,d„^ /f ) - ^ (iCu.^ [n,d„^ /f + ]- (n,a„^. In x^? + hi,) JCu,, [ff^]) • (C3) 

1=1 ^ ^ 

In the elastic case, njdn.ff^ = 6jif^^^ and the linear operator JCulfi^^] explicitly written as 

/C,,.,[/f ^ 2B,n,xf,^af,V,fl"\v). (C4) 

With this result, Eq. (|C3p becomes 



V,f'l'\ (C5) 
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Comparison with the results derived by Lopez de Haro et al. [ts'I allow us to identify liji to be defined through the 
relation 



1 / d^i 
T I an. 



1 



(C6) 



where /ij is the chemical potential of species i. 
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